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RESEARCH OVERVIEW 


Flaws exist in aircraft structures due to manufacturing operations and material defects. Under 
variable amplitude cyclic loading, these flaws grow as part-through cracks reducing the residual 
strength of structural components. To meet damage tolerant design requirements, accurate flaw 
growth predictions are needed which account for continual changes in crack shape as well as 
crack growth retardation and acceleration. Predicting the growth of part-through cracks under 
cyclic loading using an innovative and computationally efficient model is the focus of the 
research summarized in this report. 



FIG. 1 Surface crack or part-through crack. 

For through-cracks, both yield zone models and models based on plasticity-induced closure have 
been used to simulate crack growth under variable amplitude loading. Closure based models 
such as the modified strip-yield model incorporated within FASTRAN enable the prediction of 
crack growth behavior for a wide variety of load spectra without excessive empirical adjustment. 
Closure based models also allow consideration of small crack effects as recent studies suggest 
that closure plays a key role in the growth of small flaws. However, such models do not 
currently exist for the commonly encountered part-through crack shown in Fig. 1 . 

In this research effort, a slice synthesis methodology was developed and used to construct a 
modified strip-yield model for the part-through semi-elliptical surface flaw, enabling prediction 
of plasticity-induced closure along the crack front and subsequent fatigue crack growth under 
constant amplitude and variable amplitude loading. While modeling the plasticity-induced 
closure in a part-through flaw may be performed using three dimensional elastic-plastic finite 
element analysis, this type of effort is impractical from an engineering perspective. A modified 
strip-yield model similar to that used in FASTRAN for part-through flaws is a much needed 
engineering design tool, particularly when computational resources are limited. The model 
requires the following inputs: 
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1 . material modulus, Poisson’s ratio, and flow stress 

2. crack growth rate material characterizations da/dN =j{AK e fj) and dc/dN = g(AK e ff) 

3 . constraint factors at the deepest point of penetration a\ and the free surface «r 

4. plate width 2 W and thickness t 

5 . distribution of applied stress in the uncracked plate 

6 . initial flaw size specification a and c 

7. load history 

With these inputs, the developed strip-yield model gives: 

1 . crack opening stresses at the deepest point of penetration and the free surface as a function 
of the number of applied cycles N 

2. flaw depth a and length c as a function of the number of applied cycles N 

The slice synthesis methodology used relies on approximate weight functions for two- 
dimensional slices which in turn approximate the three-dimensional surface crack geometry as 
illustrated in Fig. 2. The weight functions associated with the slices shown are approximated as 
an interpolation between a fixed condition with no displacement where the spring stiffness k -> 
oo and a free condition with k-> 0. To achieve the required interpolation, weight functions for 
both the free m 0 and fixed conditions m x were needed, with the interpolation written as 

m = m o0 + rj{m 0 - ) 

where 77 is a transition factor which varies form zero to unity. An expression of this type is 
required for both the a and c slices. 
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The transition factor for both the a and c slices was assumed to be a function of the normalized 
crack depth a/t, the aspect ratio a/c, and the normalized crack length c/W. Specific values of the 
transition factors were determined for prescribed a/t, a/c, and c/W by minimizing the following 
function using a downhill simplex method with simulated annealing 

fi/la ’ r Ic) = i^ANR ~ K A ) + BNR ~ ^ B ) 

where K ANR and K B nr are the stress intensity factors at points A and B respectively on the crack 
front as computed using the Newman-Raju equations, while K A and K B are stress intensity factors 
as computed using the slice synthesis methodology with K A = K A ( ij a ) and K B — K B (r/ c ). The 
transition factors were tabulated for numerous a/T, a/c, and c/W, and a three-dimensional linear 
interpolation was used to extract the transition factor for a/t, a/c, and c/W values of interest. 

The strip-yield model in its current state is limited to a/t < 0.4. The required weight function 
transition factors could not be successfully obtained for deeper cracks. 

When using the strip-yield model, simplified variable amplitude loading may be considered in 
which the loading is defined as a sequence of constant amplitude blocks of loading of prescribed 
length N, maximum stress, and load ratio R. Note that this variable amplitude loading is much 
less general than that allowed in FASTRAN. 

Details and numerous results obtained from the strip-yield model can be found in references (1), 
(2), (3), (6), and (8). 

SUMMARY OF RESULTS 

Constraint Characterization 

When using modified strip-yield models, it is assumed that the stress within the plastic zone is a 
constant. The stress <y yy in the plastic zone will actually be nonuniform due to constraint effects 
as well as strain hardening. A constraint factor is used to approximate this nonuniform stress. 
The constraint factor or describes the average stress to flow stress ratio in the plastic zone. At the 
deepest point of flaw penetration, a plane strain condition exists and a large constraint factor will 
result. At the free surface, a plane stress condition exists and a smaller constraint factor will be 
observed. 

Elastic-plastic finite element analyses of surface-cracked plates were performed using the 
commercial finite element code ANSYS. Various crack geometries were analyzed under tension 
and bending loads. A constraint factor dependent on the location along the perimeter of the 
surface crack was computed. This constraint factor was the averaged normal stress to flow stress 
ratio acting on a line originating on and perpendicular to the crack front at a prescribed location 
and terminating at the perimeter of the plastic zone on the crack plane. 

The finite-element meshes were generated using a FORTRAN code mesh3d scpcell written by J. 
Faleskog (Royal Institute of Technology, Stockholm, Sweden). The size of the models used 
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ranged from about 33,500-51,500 nodes and 30,000-46,400 eight-noded hexahedral elements. A 
typical model is shown in Fig. 3. It should be noted that roughly two years of experience with 
this mesh generation code has revealed that the software has difficulty when large degrees of 
refinement are attempted, with poor element aspect ratios resulting in transition regions. 



1 c 1 

FIG. 3. Typical finite-element model of surface-crack configuration 

The constraint factor a is the average of the normal opening mode stress, <r }y , acting along a 
“ray” perpendicular to the crack front of the yielded elements. The constraint factor is defined as 



where A„ is the projected area of the yielded element n on the crack plane, (cry/a 0 )„ is the 
normalized centroidal stress perpendicular to the crack plane for element n, and A r (<f>) is the total 
projected area for all yielded elements along a given ray, Yielded elements were defined by 
a non-zero equivalent plastic strain. Each ray of elements were defined by the parametric angle, 

<f>, along the crack front. This angle was calculated at the midpoint of the two nodes of each 
element immediately adjacent to the crack front. 

Fig. 4 show some typical constraint variations along the crack front for different crack 
configurations subjected to remote bending and tension. The constraint is near unity (plane 
stress) at the free surface, as expected, and rises at the maximum depth location ( <j>= 90°). 

Details may be found in references (4) and (7). 
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Parametric angle, <f>, deg. 


FIG. 4. Typical constraint factor variation. 

Verification of Surface Flaw Modified Strip-Yield Model 

The strip-yield model for the surface flaw was partially verified by comparing monotonic plastic 
zone sizes and crack surface displacements with those computed from detailed elastic-plastic 3D 
finite element results as seen in Fig. 5. Crack surface displacements from the model were shown 
to compare well with results from three-dimensional elastic-plastic finite element analyses. 

When looking at the plastic zone on the crack plane, the plastic zone was observed to be largest 
beneath the free surface as shown in Fig. 6. However, viewing the crack front plastic zone on the 
crack plane alone results in an unrealistic perspective of crack front plasticity. A three- 
dimensional view of the part-through flaw plastic zone along the crack front revealed yielding on 
the crack plane is highly restricted when compared to regions off of the crack plane as shown in 
Fig. 7. For more detail see reference (8). 



0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1 0 

x/c z/a 

FIG. 5. Typical crack surface displacement comparison. 
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FIG. 6. Typical plastic zone variation along crack front in the crack plane. 



FIG. 7. Typical plastic zone near the free surface. 

Characterization of Fatigue Crack Growth Behavior for AA 7075-T6 

To verily the strip-yield model, experimental data describing the growth of surface cracks under 
cyclic loading was needed. 7075-T6 aluminum plate with a thickness of 0.50 inches was 
purchased for testing purposes. From this plate, surface crack fatigue crack growth specimens 
were machined as shown in Fig. 8. Multiple fatigue tests were conducted under constant 
amplitude loading with a stress ratio R = 0.1, which inherently included the effects of fatigue 
crack closure, and R = 0.7, which was assumed closure free. Crack length measurements on the 
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free surface were made with a traveling microscope. To allow crack depth measurements, 
marker bands were used. 

The aspect ratio evolutions for the two stress ratios tested were compared to determine whether 
crack closure affects aspect ratio or crack shape evolution. It was found that closure impacts this 
evolution in a significant manner as shown in Fig. 9. The observed crack growth rate-stress 
intensity factor range relationships for R = 0.7 at the free surface and the deepest point of 
penetration were compared to existing through-crack data, with a good correlation being 
observed as shown in Fig. 10. 

In addition, the fatigue crack growth data collected from the multiple specimens was used to 
perform uncertainty analyses on the crack lengths, crack depths, and the aspect ratio. This was 
done to determine if crack growth models are best validated using the aspect ratio as the 
measured quantity or the individual measured crack length and crack depth. The measured 
aspect ratio was found to exhibit a smaller percentage of experimental uncertainty, and thus 
should be used for model validations. Details can be found in reference (9). 

The idea of a single intrinsic curve such as that shown in Fig. 10 for a given material in a given 
environment is an integral part of fatigue crack growth prediction methodologies incorporating 
crack closure concepts. The existence of such a curve has come into question, as recent research 
has suggested that da/dN is a function of both AK e ff and the maximum stress intensity K max . 

During the next stage of the research, which is currently underway, fatigue crack growth rate will 
be measured in AA 7075-T6 compact tension specimens under both K max control and load 
control allowing determination of da/dN as a function of AK eff in both the Paris regime and the 
threshold regime. In the threshold and near threshold regime, load shedding with multiple fixed 
Kmax levels will be used to ascertain effects of K max on the da/dN - AK e g curve. In the Paris 
regime, opening load measurement will be improved using a newly developed smoothing 
operation as discussed in reference (5). Given that the fatigue crack growth process is highly 
stochastic in nature, multiple tests will be performed using each of the K max values, and a detailed 
experimental uncertainty analysis will be performed. This will allow uncertainty bands for the 
data to be established. Any observed K max effects can only be considered significant if a change 
in K max produces data outside the uncertainty bands. 



thickness = 12.7 50.80 

all dimensions are mm 


FIG. 8. Surface crack fatigue crack growth specimen. 
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FIG. 9 Aspect ratio evolution. 



FIG. 10. AA 7075-T6 fatigue crack growth rate for through-cracks and surface cracks. 

Modeling Growing Surface Flaws Using Three-Dimensional Elastic-Plastic Finite Element 
Analysis 

Crack opening stresses for surface cracks were also predicted using finite element analyses to 
allow comparison with the strip-yield model. This effort is not yet completed. Crack opening 
levels were determined for two different surface crack models using ANSYS. The first surface 
crack model was a semi-circular crack with an aspect ratio, a/c=\. The second model was a 
semi-elliptical surface crack with an aspect ratio, ale - 0.5. 
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The semi-circular model shown in Fig. 1 1 was used to determine the effects of various solution 
options. Three different solution options were studied: the loading increment, the iterative 
solver tolerance, and the use of the non-linear geometry (large deformation) option. The model 
used a relatively coarse mesh, so high stress levels ( S max / cr 0 = 0.7) were required to ensure 
adequate refinement in the plastic zone. At the maximum load, the model had 1 5 elements in the 
plastic zone at the crack deep point, and 33 elements in the plastic zone near the free surface. 

The mesh had a total of 3 1 ,392 nodes and 27,768 elements. The crack opening levels along the 
crack front were determined after 10 crack growth cycles. Four analyses were run with this mesh 
to determine the best solution options. All four analyses resulted in the same opening levels 
along the crack front (Fig. 12), so solution time can be used to determine the best options. The 
solution options and time are shown in the table below. It can easily be seen from these results 
that the best solution options in terms of time should be with a loading increment of 5% of the 
load range, a solver tolerance of IE-4, and large deformation effects off. A final analysis will be 
run with these options used together to insure there is no effect on the crack opening levels. 



Load Increment 

PCG Solver 
Tolerance 

Large 

Deformations 

Solution Time 
(hours) 

Case 1 

2.5% 

IE-8 

OFF 

71 

Case 2 

5% 

IE-8 

OFF 

48.5 

Case 3 

2.5% 

IE-4 

OFF 

57 

Case 4 

2.5% 

IE-8 

ON 

92 


A semi-elliptical model shown in Fig. 13 was used to compare finite element results with 
published experimental results. The model geometry was chosen to match a specimen tested by 
Putra and Schijve. A large model was required to ensure adequate mesh refinement in the plastic 
zone. At maximum load there were 1 0 elements in the plastic zone at the crack deep point. The 
mesh contained 59,493 nodes and 53,804 elements. Since the model was so large, the crack 
opening levels were determined after only 5 crack growth cycles. Because the crack was grown 
only 5 cycles, the opening levels at the free surface did not stabilize, but this is not of too much 
concern since opening levels near the free surface were not reported by Putra and Schijve. The 
opening levels for the crack interior as determined by FEA are much higher than those 
determined experimentally as seen in Fig. 14. This effort is currently continuing. 
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